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Abstract. We develop the analogy between self-gravitating Brownian particles and bacte- 
rial populations. In the high friction limit, the self-gravitating Brownian gas is described by 
the Smoluchowski-Poisson system. These equations can develop a self-similar collapse leading 
to a finite time singularity. Coincidentally, the Smoluchowski-Poisson system corresponds to a 
simplified version of the Keller-Segel model of bacterial populations. In this biological context, it 
describes the chemotactic aggregation of the bacterial colonies. We extend these classical models 
by introducing a small-scale regularization. In the gravitational context, we consider a gas of 
self-gravitating Brownian fermions and in the biological context we consider finite size effects. In 
that case, the collapse stops when the system feels the influence of the small-scale regularization. 
A phenomenon of "explosion" , reverse to the collapse, is also possible. 
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1. Introduction. Self-gravitating systems such as globular clusters and elliptical 
galaxies form a Hamiltonian system of particles in interaction that can be supposed 
isolated in a first approximation Since energy is conserved, the proper statistical de- 
scription of stellar systems is the microcanonical ensemble 0. The dynamical evolution 
of elliptical galaxies is governed by the Vlasov-Poisson system which corresponds to a 
collisionless regime. On the other hand, the kinetic theory of stars in globular clusters 
is based on the Landau equation (or the orbit averaged Fokker-Planck equation) which 
describes a collisional evolution. These equations conserve mass and energy. Furthermore, 
the Landau equation increases the Boltzmann entropy (H-theorem) due to stellar encoun- 
ters. These equations have been studied for a long time in the astrophysical literature and 
a relatively good physical understanding has now been achieved. In particular, globular 
clusters can experience core collapse related to the "gravothermal catastrophe" 

For systems with long-range interactions, statistical ensembles are not equivalent. 
Therefore, it is of conceptual interest to compare the microcanonical evolution of stel- 
lar systems to a canonical model. This can be achieved by considering a gas of self- 
gravitating Brownian particles submitted to a friction with an inert gas and a stochastic 
force, in addition to self-gravity 0). This system has a rigorous canonical structure. In 
the mean-field approximation, the self-gravitating Brownian gas model is described by 
the Kramers-Poisson system. In a strong friction limit, or for large times, it reduces 
to the Smoluchowski-Poisson system. These equations conserve mass and decrease the 
Boltzmann free energy. They possess a rich physical and mathematical structure and 
can lead to a situation of "isothermal collapse" 0, which is the canonical version of the 
"gravothermal catastrophe". These equations have not been considered by astrophysi- 
cists because the canonical ensemble is not the correct description of stellar systems and 
usual astrophysical bodies do not experience a friction with a gas (except dust particles 
in the solar nebula Yet, it is clear that the self-gravitating Brownian gas model is of 
considerable conceptual interest to understand the strange thermodynamics of systems 
with long-range interactions and the inequivalence of statistical ensembles. 

In addition, it turns out that the same type of equations occur in biology in relation 
with the chemotactic aggregation of bacterial populations 0. A general model of chemo- 
tactic aggregation has been proposed by Keller & Segel |H1 in the form of two coupled 
differential equations. In some approximation 9 , this model reduces to the Smoluchowski- 
Poisson system, exactly like for self-gravitating Brownian particles. Therefore, there exists 
an isomorphism between self-gravitating Brownian particles and bacterial colonies. In this 
paper, we shall develop this analogy in detail. We shall also propose a modification of 
the "standard model" by introducing a small-scale regularization. In the gravitational 
context, we shall invoke Pauli's exclusion principle and consider a gas of self-gravitating 
Brownian fermions. In the biological context, we shall heuristically account for finite size 
effects by considering a lattice model. In that case, the collapse stops when the system 
feels the small-scale regularization. An explosion phenomenon, reverse to the collapse, is 
also possible. Finally, we shall discuss the difference between elliptical and parabolic mod- 
els of bacterial populations and gravitational systems. We shall also show that vortices 
in two-dimensional turbulence exhibit features similar to stars and bacteries. 
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2. The Hamiltonian A^-stars problem Consider a system of N stars in gravita- 
tional interaction. We assume that the system is isolated so that it conserves mass and 
energy. The equations of motion can be cast in a Hamiltonian form 

dvi dH dvi dH 

dt ovi dt OTi 



where H is the Hamiltonian 



^ ^=l t<3 "^Jl 

This A^-body problem is the correct starting point in the description of globular clusters 
and elliptical galaxies Since the system is isolated, the relevant statistical ensemble 
is the microcanonical ensemble |2] . The statistical mechanics of self-gravitating systems 
enclosed within a spherical box of radius R was initiated by Antonov |10| and Lynden- 
BeU & Wood 0. They found that for A = -ER/GPP > 0.335, corresponding to smaU 
energies, there is no maximum entropy state so that the system must collapse (see Fig. 
This is the so-called "gravothermal catastrophe". 

In the case of elliptical galaxies, encounters between stars are completely negligible 
for the timescales of interest and the dynamics is described by the Vlasov-Poisson system 

^ + v — -FF- — = f3) 
dt dr dv ' 

A$ = AnG J fd^v, (4) 

where /(r, v,t) is the distribution function, F(r,t) = — V$ is the gravitational force and 
<i>(r,i) is the gravitational potential. By contrast, for globular clusters, stellar encounters 
must be taken into account. In that case, the dynamics is described by the Landau 
equation 

Tt^^'-^^^-d^-d^^j'^^'^ t^^9^"^9</' 
where K^^'^ is the tensor 



A 



K^- ^ (6) 

u \ J 

where /i = /(r, vi, i), u = vi - v is the relative velocity and A = 27rG^mln(L 

max I ■^min) 

is a constant {Lmax and Lmin are appropriate lengthscales). This equation conserves 
mass M = J pd^r and energy E = J f^d^rd'^v + ^ J p^d^r and satisfies a i7-theorem 
[S > 0) for the Boltzmann entropy 

= - //In/d^rdV (7) 



Therefore, due to the development of encounters between stars, the system is expected 
to relax towards a statistical equilibrium state described by the Maxwell-Boltzmann 
distribution 

/ = Ae-'^('#+*), (8) 
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Figure 1: Caloric curve for self-gravitating systems exhibiting a "gravothermal catas- 
trophe" in the microcanonical ensemble and an "isothermal collapse" in the canonical 
ensemble. 

which maximizes the Boltzmann entropy at fixed mass and energy (microcanonical de- 
scription). In fact, this is the case only if the energy is sufficiently high. Below the Antonov 
critical energy, there is no maximum entropy state and the system undergoes core col- 
lapse. This is a manifestation of the "gravothermal catastrophe" . Dynamical models show 
that the collapse is self-similar and leads to a finite time singularity (the central density 
becomes infinite in a finite time). By solving the orbit averaged Fokker- Planck equation 
numerically, Cohn 11 finds that the density profile behaves as p ~ r"" with a — 2.23 
at large distances. Alternatively, Lancellotti & Kiessling V2, consider the full Landau- 
Poisson system and argue that a = 3. 

3. The self-gravitating Brownian gas The Hamiltonian equations Q-© describe 
an isolated system for which energy is conserved. It can be of interest to study in parallel 
a model that is stochastically forced by an external medium. We thus introduce a system 
of N Brownian particles in gravitational interaction described by the Langevin equations 



where — ^v^ is a friction force and Ri(i) is a white noise satisfying (Ri(f)) — and 
{Ra,i{t)Rb,j{t')) — SijSabS{t — t'), where a, 6 = 1,2,3 refer to the coordinates of space 
and i,j = 1, ...,iV to the particles. The particles interact via the gravitational potential 
?7(ri, rjv) = 'Ylii<j'^^i~'^i) where it(ri — rj) = —G/|ri — rj|. We define the inverse tem- 
perature (i — \/T through the Einstein relation ^ = D(3. The self-gravitating Brownian 
gas model has a rigorous canonical structure (see Appendices A and B) where the temper- 




(9) 
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ature T measures the strength of the stochastic force. The stochastic process Q defines 
a model of gravitational dynamics which extends the classical Einstein-Smoluchowski 
Brownian model to the case of stochastic particles in interaction. In this context, 
the friction is due to the presence of an inert gas and the stochastic force is due to clas- 
sical Brownian motion, turbulence or any other stochastic effect. This model can also be 
viewed as a generalization of the Chandrasekhar T? stochastic model which describes 
the evolution of a single test particle in a stellar cluster at statistical equilibrium (ther- 
mal bath approximation). In that context, the diffusion and the friction model stellar 
encounters. More generally, the friction and the noise can mimick the overall influence of 
an external medium (not represented) with which the particles interact. 

Starting from the A'^-body Fokker-Planck equation and implementing a mean-field 
approximation which is valid in a proper thermodynamic limit N — > +00 with rj = 
(3GM/R fixed, we show in Appendix B that the distribution function /(r,v, i) satisfies 
the Kramers-Poisson system 



df df ^ df d , ^ 

— + v ■ — +F ■ — ^ i D 

dt dr dv dv 



If+Pf^ 
dv 



(10) 



A$ = AnG J fcPv. (11) 

The Kramers-Poisson system decreases the Boltzmann free energy 

F[f] = E -TS = J /yd^rd^v+i J p<S>d^r + T J flnfd^rd^v. (12) 

Therefore, the system is expected to relax towards a statistical equilibrium state described 
by the Maxwell-Boltzmann distribution 

/ = Ae~'5('#+*), (13) 

which minimizes the Boltzmann free energy at fixed mass and temperature (canonical 
description). In fact, this is the case only if the temperature is sufficiently high. Below 
a critical temperature the free energy has no minimum and the system undergoes 
an "isothermal collapse" (see Fig. We note that the equilibrium states of stellar 
systems (e.g., globular clusters) and self-gravitating Brownian particles are both given 
by the isothermal distribution ((SJ or H13() . However, the stability limits are different in 
each ensemble because the caloric curve (3{E) presents turning points (see Fig. The 
stability of stellar systems and self-gravitating Brownian particles differs in the region of 
ensemble inequivalence where the specific heat C = dE/dT is negative 0. 

In the high friction limit ^ +00, or equivalently for large times t ^ we can 
neglect the inertia of the particles so that the Langevin equations © take the form 

e^ = -V,f/(ri,...,rjv) + ViDR,(t). (14) 

Again using a mean-field approximation (see Appendix B), we can show that the density 
p{v, t) satisfies the Smoluchowski-Poisson system 

dt 



i(TVp + pV$) 



(15) 
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T=0.2 
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Figure 2: Isothermal collapse in I? = 3. For rj = f3GM/R > 2.52, the system undergoes a 
finite time singularity leading to a p ~ density profile A Dirac peak is formed 

in the post-collapse regime |21| . 



A$ = AnGp. (16) 

The Smoluchowski equation (|15|) can also be deduced from the Kramers equation (|HJ|I 
by using a method of moments 15 or a Chapman-Enskog expansion We can show 
that it decreases the Boltzmann free energy 

F[p] J p\np (fr+^ J p<i> (fr, (17) 

obtained from Eq. (|12|l by using the fact that the distribution function is close to the 
Maxwellian distribution 

/=(^) P(r,t)e-''*, (18) 

in the high friction limit |15[ I16j. The equilibrium states of the Smoluchowski-Poisson 
system (|15ll - (|16|) are given by the Boltzmann distribution 

p = A'e-^^. (19) 

In the gravitational context, the Kramers-Poisson system and the Smoluchowski- 
Poisson system have been introduced by Wolansky J7|, Chavanis, Sommeria & Robert 
[T^ and Chavanis, Rosier & Sire @] with different motivations. A physical description of 
the spherically symmetric solutions of the Smoluchowski-Poisson system has been given in 

1191 1^ 1211 122j in various dimensions of space (including the critical dimension D = 2) 
and for both pre-coUapse and post-collapse regimes. For sufficiently large mass M or 
sufficiently low temperature T, the Smoluchowski-Poisson system displays a self-similar 
collapse leading to a finite time singularity in _D > 2 (see Fig. O. The density profile 
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Figure 3: Isothermal collapse in D = 2. For 77 = PGM > 4, the system develops a Dirac 
peak surrounded by a halo behavmg as p ^ r^". For t — + tcoii, a = 2 but the convergence 
to this asymptotic value is extremely slow so that an effective exponent a ~ 1.3 is 
observed for the times achieved in our simulations |19|. 



behaves as p ~ with a = 2 at the collapse time Then, the evolution continues 

and a Dirac peak ("black hole") is finally formed in the post-collapse regime This is 
consistent with predictions of statistical mechanics in the canonical ensemble |23l |2| ■ In 
D = 2, the evolution is more complex 19 and creates a Dirac peak containing a fraction 
T/Tc of the total mass, surrounded by a p ^ r^" halo with an effective scaling exponent 
converging very slowly to a = 2 (see Fig.l^J- A review of these results is given in 2^,. 



4. Self-gravitating Brownian fermions As discussed previously, self-gravitating 
classical particles have the tendency to develop finite time singularities. In an attempt to 
regularize the problem at high densities, and avoid unphysical infinities, we can invoke 
quantum mechanics and use Pauli's exclusion principle. Thus, we shall consider a gas of 
self-gravitating Brownian fermions as a regularized model of gravitational dynamics. The 
equilibrium states of self-gravitating fermions and the description of phase transitions 
in the self-gravitating Fermi gas (in both microcanonical and canonical ensembles) have 
been investigated by Chavanis 25 , 26 . The system of self-gravitating Brownian fermions 
|27j can be used as a simple model to study these phase transitions dynamically in the 
canonical ensemble. 

A generalization of the Kramers equation 1)10(1 taking into account the Pauli exclusion 
principle is given by |TB] : 

This equation respects the constraint f < rjo = /K^ {h is Planck constant) at all times 
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2 

a=-er/gm' 



Figure 4: Hysteretic cycle in the canonical ensemble for self-gravitating fermions. The 
system undergoes a "collapse" at Tc and an "explosion" at T*. The temperatures of 
collapse and explosion differ due to the existence of long-lived metastable states (local 
minima of free energy) \2^ . The circles correspond to the results of numerical simulations 
of the fermionic Smoluchowski-Poisson system |27j . 



and decreases the Fermi-Dirac free energy 

^ " / Y^''"^'^ +IJ ^'^^'^ + T y {/ In / + (770 - /) ln(??o - f)}d^rd^v. (21) 



The equilibrium states are given by the Fcrmi-Dirac distribution function 

/ 



(22) 



1 + Ae'3(^+*) 

which minimizes the free energy (|21|l at fixed mass and temperature. 

Now, considering the high friction limit |15[ , we can derive a generalized Smolu 
chowski equation of the form 



dp 
dt 



(Vp + pV$) 



where p{p) is the local equation of state of the Fermi gas 



^3/2 



and 



In{t) 



87rV2?7o 



-dx 



(23) 



(24) 



(25) 



1 + ie^ 

is the Fermi integral. The fermionic Smoluchowski-Poisson system decreases the free 
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Figure 5: Collapse of self-gravitating Brownian fermions for T < Tc. The system finally 
forms a "fermion ball" (similar to a white dwarf star) surrounded by an isothermal halo 
(like a vapor). This is a global minimum of free energy in the canonical ensemble. 



F[p]^ I pl '^dp'd^r+i /p$d^r, (26) 



energy 



" 2 

obtained from Eq. (|21|l by using the fact that f{r,v,t) is close to the Fermi distribu- 
tion in the high friction limit ^H!- Furthermore, the equilibrium states are obtained by 
substituting the relation 

P- ^3/2 h^^'' )' (27) 

in the Poisson equation p6(l and solving for the gravitational potential $ |25j . 

The fermionic Smoluchowski-Poisson system has been studied mathematically by Biler 
et al. j28j . Numerical simulations have been conducted in parallel in |27j . In particular, 
these simulations reveal an interesting hysteretic cycle (see Fig. 0J) discussed by Chavanis 
& Rieutord [221 . Below a critical temperature Tc, coinciding with the Jeans instability 
criterion [HIH], the system collapses under its own gravity as in Fig. El However, for self- 
gravitating fermions, the collapse stops when the core becomes degenerate (in the sense 
of the Fermi-Dirac statistics). The resulting condensed object has the same structure as 
a cold white dwarf star (or a "fermion ball") in which gravity is balanced by quantum 
pressure (see Fig.lSJ. If now temperature is increased, the system remains in the condensed 
phase until another temperature T* > Tc (see Fig. 0)) at which it undergoes an explosion 
(reverse to the collapse) and returns to the gaseous phase (see Fig. E)). 
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Figure 6: Explosion of self-gravitating Brownian fermions for T > T^,. When diffusion 
prevails over gravity, the condensed object expands and the system returns to the gaseous 
phase. 



5. Chemotactic aggregation We shall now point out some analogies between a 
gas of self-gravitating Brownian particles and the chemotaxis of bacterial populations [7] . 
The name chemotaxis refers to the motion of organisms (amoeba) induced by chemical 
signals (acrasin). In some cases, the biological organisms secrete a substance that has 
an attractive effect on the organisms themselves. Therefore, in addition to their diffusive 
motion, they move systematically along the gradient of concentration of the chemical 
they secrete (chemotactic flux). When attraction prevails over diffusion, the chemotaxis 
can trigger a self-accelerating process until a point at which aggregation takes place. This 
is the case for the slime mold Dictyostelium Discoideum and for the bacteria Escherichia 
coli. 

A model of slime mold aggregation has been introduced by Keller & Segel [Sj in the 
form of two PDE's: 

|^=V-(i?2Vp)-V-(i^iVc), (28) 
dc 

— = -k{c)c + f{c)p + D,Ac. (29) 
ot 

In these equations p{r, t) is the concentration of amoebae and c(r, t) is the concentration 
of acrasin. Acrasin is produced by the amoebae at a rate /(c). It can also be degraded 
at a rate k{c). Acrasin diffuse according to Pick's law with a diffusion coefficient D^. 
Amoebae concentration changes as a result of an oriented chemotactic motion in the 
direction of a positive gradient of acrasin and a random motion analogous to diffusion. 
In Eq. (|28|l . D2{p, c) is the diffusion coefficient of the amoebae and Di{p, c) is a measure 



11 



of the strength of the influence of the acrasin gradient on the flow of amoebae. This 
chemotactic drift is the fundamental process in the problem. 

A first simplification of the Keller-Segel model is provided by the system of equations 

|^=i?Ap-xV-(pVc), (30) 
dc 

— =D'Ac + ap-bc, (31) 

where the parameters are positive constants. An additional simplification, introduced by 
Jager & Luckhaus consists in ignoring the temporal derivative in Eq. H31() . This is 
valid in the case where the diffusion coefficient D' is large. Taking also 6 = 0, we obtain 

^^DAp-xV -ipVc), (32) 

Ac = -A/9, (33) 

where A = a/D'. Clearly, these equations are isomorphic to the Smoluchowski-Poisson 
system (|15|) - (|16|l describing self-gravitating Brownian particles in a high friction limit. 
In particular, the chemotactic flux plays the same role as the gravitational drift in the 
overdamped limit of the Brownian model. We have the correspondance $ ^ — ^^c, 13 ^ 
^^Qjy and ^ ^ . Through this analogy, we can develop an effective thermodynamical 
formalism to investigate the chemotactic problem ^11301 . In particular, Eq. (|32|) is similar 
to a Fokker-Planck equation and its stationary solutions are similar to the Boltzmann 
distribution 

p = Ae*^ (34) 
which maximizes a Lyapunov functional similar to the Boltzmann free energy. 

6. A generalized chemotactic model The Keller-Segel model ignores clumping 
and sticking effects. However, at the late stages of the blow-up, when the density of 
amoebae has reached high values, finite size effects and stickiness must clearly be taken 
into account. As a first step, we propose to replace the classical equation (|32|) by an 
equation of the form 

|^ = i?Ap-xV-(p(l-p/ao)Vc), (35) 

which enforces a limitation p < ag on the maximum concentration of bacteria in physical 
space 31 . This equation increases the Lyapunov functional 

J[p] = - J[p\np+ (no - p) Hao - p)]d^r + ^ J d^^. (36) 

In the thermodynamical analogy mentioned above, this functional can be interpreted 
as a free energy J = S — f3E associated with a Fermi-Dirac entropy in physical space 
|15[ I3U| . This form of entropy can be obtained by introducing a lattice model preventing 
two particles to be on the same site. The lattice creates an exclusion principe in physical 
space similar to the Pauli exclusion principle in phase space. Then, S[p\ can be obtained 
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by a standard combinatorial analysis respecting this exclusion principle. The equilibrium 
states of Eq. H35|l are given by a Fermi-like distribution in physical space 

P = — (37) 

1 + Ae o 

which maximizes the effective free energy (|36|l at fixed mass. 

It is clear on qualitative grounds that the system of equations H35 |l - H33(> will display 
exacty the same phenomena as the system of self-gravitating Brownian fermions (at least 
in dimension D — i) (31| . In particular, by tuning the mass of the system, or more 
generally the dimcnsionless parameter -q = a-kDR ' '^^'^ describe an hysteretic cycle 
similar to the one depicted in Fig. 21 Above a critical mass, blow-up occurs until finite 
size effects come into play and arrest the collapse. Then, by slowly decreasing the mass of 
the aggregate, we reach a critical point 77, at which an explosion sets in. This occurs when 
diffusion prevails over chemotactic drift. This hysteretic cycle has never been reported in 
the chemotactic literature and we can wonder whether it could be observed in biological 
experiments. In any case, the constraint p < ctq implied by Eq. H35|l regularizes the 
problem and prevents unphysical infinities. 

On the other hand, as pointed out by Keller & Segel |H], the diffusion coefficient of 
amoebae can depend on the density p, leading to a situation of anomalous diffusion. 
For example, the case where the diffusion coefficient is a power law of the density has 
been investigated in [201 ■ Moreover, the relation between the concentration of amoebae 
and acrasin may be more complex that simply given by the Poisson equation 1)33(1 . For 
example, taking 6 7^ in the original model, we obtain a relation of the form 

Ac-/c^c=-Ap, (38) 

where A:| = b/D' . The second term is similar to the Debye shielding in plasma physics. 
These remarks motivate us to consider a larger class of drift-diffusion equations of the 
form 

dp 



{^[Vp+^V j u{v~v')p{v',t)d^v' 



(39) 



This generalized class of Fokker-Planck equations has been introduced by Chavanis 
E^lEllSOl- They include an arbitrary diffusion coefficient D{r,t), an arbitrary convex 
function C{p) and an arbitrary binary potential of interaction w(r — r'). Equation H39II 
can therefore provide a generalized model of chemotactic aggregation taking into account 
anomalous diffusion, stickening effects and shielding effects of various forms. This equation 
increases a Lyapunov functional 

J[p] = - / Cip)d^r - i/3 / p(r, t)u{r - r')p{r', t)d^rd^r', (40) 



2 

which plays the role of a generalized free energy in an effective thermodynamical formal- 
ism. Finally, the stationary solutions of this equation are determined by the intcgrodif- 
ferential equation 

C'{p) ^-P J uir - r')p(r')d^r - a. (41) 
In the limit of short range interactions, the non-local drift-diffusion equation H39|l reduces 
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to a form of Cahn-Hilliard equation 

|^=V.|xV(Ap-^^'(p))|, (42) 

where V{p) = -(2/h(i)C{p) - {a/h)p^, x - /3|c^ ^i*^ " ^ S and b = 

/ u(|x|)a;^(i^x POIEI- The Cahn-HiUiard equation has been extensively studied in the 
theory of phase ordering kinetics. In the chemotactic model, the short-range interaction 
limit is reached when fc 3> 1, i.e. in a regime of high degradation rate. The Cahn-Hilliard 
equation is known to develop "domain walls" and other morphological structures. It 
would be interesting to see whether such solutions can be constructed and observed in 
the context of bacterial aggregation. 

The drift-diffusion equation H39|) does not take into account memory effects. However, 
if we come back to the original Keller-Segel model, the concentration of acrasin c is related 
to the concentration of amoeba p by an equation of the form 

3c 

— = D'Ac + ap- be, (43) 
at 

which involves a time derivative. Therefore, a more general model of chemotaxis is rep- 
resented by the non-Markovian equation 



C"{p) 

taking into account delay effects 



Vp+tJ-tV / f\{r^r',t-t')p{r',t')d''r'dt' 







(44) 



7. Hyperbolic models Recent experiments of in vitro formation of blood vessels 
show that cells randomly spread on a gel matrix autonomously organize to form a con- 
nected vascular network that is interpreted as the beginning of a vasculature. This phe- 
nomenon is responsible of angiogenesis, a major actor for the growth of tumors. These 
networks cannot be explained by the above parabolic models that lead to pointwise blow- 
up. However, they can be recovered by certain hyperbolic models. In particular, Gamba 
et al. [HS! describe the evolution of the cell population by the equations 

dp dvi 

^ V • (pu) = 0, — u • Vu = -V/i(p) -t- p\7c, (45) 

— = D'Ac + ap- be, (46) 
ot 

and show, through extensive numerical simulations, that these equations develop network 
patterns that are in good agreement with experimental results. Additional numerical 
simulations are performed in (see Fig.O. 

It can be noted that these equations are similar to the Euler-Jeans equations 

1^ + V • (pu) = 0, ^ + u • Vu = ~Vh{p) - V$, (47) 



A$ = inGp, 



(48) 
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Figure 7: Formation of a network in the hyperbolic model H45() - (|46ll of chemotactic ag- 
gregation pi]. We kindly acknowledge F. Filbet for having provided this figure. 

considered in the astrophysical literature (in that context, h{p) is the enthalpy). These 
equations also describe the emergence of network-like patterns in the large-scale distri- 
bution of masses in the universe (in that case we need to account for the expansion of 
the universe). This is another aspect of the beautiful analogy between astrophysics and 
biology. 

The connexion between parabolic and hyperbolic models can be made explicit by 
introducing a friction force — in the momentum equation. This friction term can take 
into account the influence of the gel matrix on the motion of cells in biology or the friction 
with a gas in the case of dust particles moving through the solar nebula ^ . In the high 
friction limit, or for sufficiently large times t ^ the inertia of the particles can be 
neglected and the parabolic models 1)1511 and H32|l are recovered. 



8. The Hamiltonian iV- vortex problem We finally conclude this discussion by 
showing that two-dimensional vortices in fluid mechanics also share some analogies with 
self- gravitating systems and bacterial populations. Consider a system of N point vortices 
in 2D hydrodynamics. The equations of motion can be cast in a Hamiltonian form 

dxi dH dy dH 



dt diji ' dt dxi 



where H is the Hamiltonian 



i/ = -^^7^1n|r.-r 



(49) 



(50) 



and 7 the circulation of a point vortex. Since the total energy is conserved, the correct 
statistical description is the microcanonical ensemble. The statistical mechanics of point 
vortices was initiated by Onsager 35 . He showed the existence of negative temperature 
states /3 < at which point vortices cluster in macrovortices. This property is related to 
the emergence of large-scale vortices in atmospheric flows like, e.g., Jupiter's great red 
spot. 
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The statistical equilibrium state of a system of point vortices has been determined by 
Joyce & Mongomery in a meanfield approximation. It is obtained by maximizing the 
Boltzmann entropy 



S = - j ujlnuod^v, (51) 

at fixed circulation F = J ijJ(Py and energy E — ^ J uj'ilxpY, where is the vorticity and 
tp the streamfunction. This leads to the Boltzmann distribution 

uj^Ae-l^'^', (52) 

where the stream-function ?/; is determined self-consistently by substituting the relation 
in the Poisson equation 

uj = -Ai/), (53) 

and solving the resulting Boltzmann-Poisson system. We note that the equilibrium dis- 
tribution of point vortices at negative temperatures (|52|) (|53|l is similar to the equilibrium 
distribution of stars (O (01) and bacteries H34() H33() . Therefore, there is a deep analogy be- 
tween two-dimensional vortices, stellar systems and biological clusters as first emphasized 
in EH ESI. 

A kinetic theory of point vortices has been developed by Chavanis , using methods 
inspired from stellar dynamics and plasma physics. Starting from the A'^-body Liouville 
equation, using projection operator technics and considering axisymmetric flows, he ob- 
tained a kinetic equation of the form 

dijj 7 d 



dt 4r dr 



ridri6{Vl ~ In 



1 



1 duj 1 doji 1 

-UJl- CJ— \, 

r or ri ori J 



(54) 

where fl — il{r,t), fli — n{ri,t) and r> (resp. r<) is the biggest (resp. smallest) of r 
and ri. This equation governs the evolution of the average vorticity w(r,t) = "/p(r,t) 
which is proportional to the local density p{r, t) of point vortices. The angular velocity 
r2(r, <) = {Ve){r,t)/r is related to the vorticity by 

co = --^{nr'). (55) 
r Or 

More general kinetic equations, accounting for memory effects, have also been derived 
in but they are difficult to study. Equation l|54|) is the counterpart of the Landau 
equation Q in gravitational dynamics. The evolution of the vortex cloud is due to a 
condition of resonance encapsulated in the (5-function. This resonance takes into account 
distant collisions between point vortices and is manifest only if the profile of angular 
velocity is non-monotonic (otherwise the collision term vanishes). The kinetic equation 
(j54|l conserves all the constraints of the A^- vortex system (circulation, energy, impulse and 
angular momentum) and increases the Boltzmann entropy (|51|l . However, as discussed 
in 38 , this kinetic equation does not relax towards the thermal equilibrium state H52(l 
predicted by the statistical mechanics of Joyce & Montgomery 36 . The reason is not 
well- understood. This may be a signature of the lack of ergodicity in the point vortex 
gas. This may also correspond to the break-up of the assumptions leading to Eq. H54|l . 
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Non-trivial three body correlations may be necessary to produce further evolution |88j . 
However, these correlations will be manifest on longer timescales and the system can 
remain frozen in a sort of metastable state for a very long time. 

9. Brownian vortices As in the case of self-gravitating systems, it is of concep- 
tual interest to compare the microcanonical evolution of a Hamiltonian system of point 
vortices to a canonical model. We thus introduce formally a system of N "Brownian 
vortices" (the counterpart of self-gravitating Brownian particles) interacting through a 
logarithmic potential in two dimensions. By definition, this system is described by the 
Langevin equations 

— i = -z X V,{/(ri, ...,rjv) ~ eV,C/(ri, r^) + V2DR,{t), (56) 
at 

where — ^ViC/(ri, tat) is a drift and Ri(t) is a white noise. The potential of interaction 
is given by C/(ri, ...,rjv) = Y.i<j H^i^'^j) where u{ri-rj) = -^ln|ri-rj|. Finally, we 
define the inverse temperature (3 ~ 1/T through the Einstein relation ^ = Df3. Since the 
temperature T is fixed, the Brownian vortex model has a rigorous canonical structure. 
This model can be seen as a generalization of the stochastic model of Chavanis [3^1 EHl 
which describes the evolution of a single test vortex in a bath of field vortices at statistical 
equilibrium (this is the counterpart of the Chandrasekhar ^1] stochastic model in stellar 
dynamics). In that case, the drift and the diffusion arise from discrete interactions with 
the field vortices which play the role of a thermal bath. Equation H56|l without the first 
term can also provide a stochastic model of chemotaxis. 

Starting from the A^-body Fokker-Planck equation, implementing a mean-field ap- 
proximation and following the lines of Appendix B, we find that the average vorticity 
uj{r, t) = 'yp{r, t) satisfies a Fokker-Planck equation of the form 

^-hu-VtJ = V- [£i(Vcj-f /37WVV')], (57) 
ot 

coupled to the Poisson equation 

Ai/) = -uj. (58) 

The drift-diffusion equation H57|l conserves circulation and increases the free energy 
J ^ S — PE (canonical description). Moreover, the Boltzmann distribution H52|l is the 
only stationary solution of this equation contrary to the kinetic equation H54I) valid for 
the Hamiltonian system. We note the analogy with the Smoluchowski-Poisson system 
(jl5|l - H16|l and (|32|l - 1)33(1 describing self-gravitating Brownian particles and bacterial pop- 
ulations. In the case of Brownian point vortices, the Fokker-Planck equation directly 
takes the form of a drift-diffusion equation while for self-gravitating Brownian particles, 
the Smoluchowski equation is obtained in a high friction limit where the motion of the 
particles is overdamped. We also emphasize the difference of behaviours between the 
Brownian vortex model 1)56(1 - ((57(1 and the Hamiltonian vortex model ((49(1 - ((54(1 . In the 
canonical ensemble in _D = 2, there exists a critical inverse temperature j3c = below 
which a system of Brownian vortices collapses 19 . By contrast, in the microcanononical 
ensemble, there is no collapse uvD — 2. This is related to the fact that, in two dimensions, 
there exists a maximum entropy state for all values of energies while the free energy has a 
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minimum only for (3 > Pc (the states with (3 < (3c are not accessible in the microcanonical 
ensemble) j37j . Therefore, a naive canonical description of point vortices in two dimen- 
sions does not describe the true Hamiltonian system correctly. This is a manifestation of 
ensemble inequivalence for systems with long-range interactions, such as point vortices. 

10. Violent relaxation and metaequilibrium states In the thcrmodynamical 
limit N — !■ -\-(yo with 77 = PN^j"^ fixed, appropriate to point vortices, the collision term 
appearing in the kinetic equation (|54(l is of order N^^. It represents therefore the first 
order correction to the Vlasov limit 03- However, in general, the "collisions" are negli- 
gible in the point vortex gas and the evolution of the vorticity/density is described for 
timescales of interest by the 2D Eulcr-Poisson system 

^-fu-Vw = 0, (59) 
at 

Aijj = -w. (60) 

This corresponds to the Vlasov regime in stellar dynamics and plasma physics. These 
equations can be derived from the Liouville equation by assuming that the TV-body dis- 
tribution function factorizes as N one-body distribution functions, which is justified in 
the limit N +00 (for fixed t) The Euler equation also describes the inviscid 

evolution of continuous vorticity flows in two dimensions (see, e.g., [37j). 

It is known that the 2D Euler-Poisson system develops a complicated mixing process 
leading to a metaequilibrium state on the coarse-grained scale. This process is called vio- 
lent relaxation. A small-scale parametrization of the 2D Euler equation can be attempted 
by using thcrmodynamical arguments j40| and kinetic theory |41| . Similar arguments can 
be developed for the Vlasov-Poisson system jJS] . In the two-levels approximation of their 
statistical theory, Robert & Sommeria j4(J| have proposed a parametrization of the form 

^ -l-u- Vz:J = V • \D{Vlj + (3{t){uJ ^ a-i)(ai ~ uj)Vij)], (61) 
at 

where ZU and ip are coupled through the Poisson equation l|60|l and (3{t) evolves with 
time so as to satisfy the conservation of energy E = 0. This equation respects the con- 
straints cr_i <itJ<(yi imposed by the Euler equation on the coarse-grained vorticity and 
converges towards the equilibrium distribution 

which maximizes the mixing entropy 

S = - J[p\np+il-p)\n{l- p)]d^r, (63) 

with lU = pai + (1 — p)cr_i, at fixed circulation and energy. An example of relaxation 
towards statistical equilibrium is shown in Fig. |Hlin relation with the formation of large- 
scale vortices in the jovian atmosphere 021 ■ We note the analogy between Eq. H61|l and the 
model (|35|l proposed for the chemotactic aggregration of bacterial populations (consider 
in particular the case cr_i = 0). One important difference is that (3 in Eq. (|61|) is not 
fixed but evolves with time so as to satisfy the conservation of energy. Hence, Eq. (|61l) 
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Figure 8: Relaxation towards statistical equilibrium in a QG model of Jupiter's great red 
spot and other jovian vortices (white ovals), using the parameterization H()l|l j42| . We 
kindly acknowledge F. Bouchet for having provided this figure. 



describes a microcanonical situation while Eq. H35|l can be associated with a canonical 
situation ^1 [30] . The well-posedness of Eq. <(HT)l has been established in 031 ■ 

In the quasi-geostrophic approximation appropriate to geophysical fluid dynamics [IS] , 
Eqs. H59|) and H61I) remain valid provided that the vorticity uj is replaced by the potential 
vorticity q which is related to ?/; by a relation of the form 



Kv). 



(64) 



In this equation, h{y) represents a topography and R is the so-called Rossby radius. 
It takes into account the deformation of the surface of the fluid. This equation should 
be compared with the similar equation (|38|) in the chemotactic problem. In the limit of 
small Rossby radius i? 0, the relaxation equation (|61|l becomes equivalent to the Cahn- 
Hilliard equation as discussed in Ref. In that case, its stationary solutions describe 
"domain walls" . These solutions precisely account for the jet structure of Jupiter's great 
red spot 

The parametrization of Robert & Sommeria j4()| is complicated to implement in the 
case of realistic initial conditions with a large number of vorticity levels. In order to 
go beyond the two-levels approximation while leaving the problem tractable, Chavanis 
[151 0ni has proposed an equation of the form 



dq 
'Ft 



Vq = V ■ ID 



Vq 



m 

C"{q) 



(65) 
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where C is a convex function. The relation between the potential vorticity and the stream- 
function is written in the general form 

ip{T,t)^ j u{v-v')q{r',t)d?v'. (66) 

The relaxation equation (|65|l provides a simplified parametrization of the 2D Euler- 
Poisson system, where the function C [q] has to be adapted to the physical context (ocean 
dynamics, jovian atmosphere,...). Alternatively, it can be used as a numerical algorithm 
to construct arbitrary nonlinearly dynamically stable stationary solutions of the 2D Eu- 
ler equation, specified by a relationship C'{q) — —f3ip — a, in an attempt to reproduce 
observed phenomena |15[ BS| . 



11. Conclusion In this paper, we have emphasized the analogy between self grav- 
itating systems, two-dimensional vortices and bacterial populations. This completes the 
discussion initiated in ^] . We have introduced or quoted different types of models 
and kinetic equations. Some models have a rigorous physical justification. Others have a 
more academic or conceptual interest. The common point between these models is the 
long-range nature of the interactions leading to non-local kinetic equations or non-local 
Fokker-Planck equations EOI ■ They correspond either to a microcanonical sit- 

uation (fixed energy) or a canonical one (fixed temperature). These equations have a 
very rich mathematical and physical structure associated with phase transitions, blow-up 
phenomena and morphogenesis. It is clear that these equations have a "fundamental" 
character and that they deserve to be studied in detail. This systematic study has been 
undertaken by the authors in H UHl EOl EH EH EH • 
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Appendix A: The canonical statistical equilibrium state 

In this Appendix, we use the equilibrium BBGKY hierarchy and a mean-field ap- 
proximation to determine the statistical equilibrium state of particles in interaction in 
the canonical ensemble (a similar approach can be developped in the microcanonical 
ensemble 1481). The canonical distribution in physical space is given by 

Pjv(ri,...,r^) = ^e-'3t^(----), (67) 

where Z{P) is the partition function 

Z{(3) = J e-'''^(""i'-'""")d^ri...d^rjv. (68) 

In practice, we are not interested by the A^-body distribution function PAr(ri, tat) but 
rather by the reduced densities such as 

Pi(r)= /Pjv(r,r2,...,rjv)d^r2...d^r^, (69) 
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P2(r,r') = yPA.(r,r',r3,...,rw)d^r3...d^rAr. (70) 
Assuming that the particles are identical, the average density in r is 

p(r) = /^<5(r,-r)\=iVPi(r). (71) 



Taking the derivative of Eq. H67|l with respect to ri, we get 
dPN , . dU 



OTi OTi 

We now integrate over r2, ...,rjv to obtain 

N 



(ri,...,rw) - -/JPat — (ri,...,rw). (72) 



dPi 
dr 



i(ri) = -/3 / n^^"^^ Pn(Ji, ...,rAr) — (n, ...,yn)- 



(73) 



For a binary potential of interaction 

U{yi,...,yn) =^u{y,-Vj), M(r)=w(-r), (74) 

the foregoing expression can be rewritten 

ap N „ N ^ 



j=2 i=2 

Since the particles are identical 

N 



dPi 
9ri 



(76) 



i(ri) = -P{N - 1) / n PN{ri, ...,rjv)^(ri - ra), 
or, equivalently, 

^(ri) = -/3(7V-l) /d^ra F2(ri,r2)|^(ri-r2). (77) 
OTi J dri 

For iV 3> 1 such that /JA^w ~ 1, which defines the proper thermodynamic limit of this 
nonextensive system, it can be shown |48| that the mean field approximation is exact, 
i.e. 

P2(ri,r2) = Pi(ri)Pi(r2) + 0(l/iV). (78) 

Thus, Eq. (|77|) becomes 

l^(ri) = -f3NP,{r,)£- J d^v^ Pi(r2)u(ri - v^). (79) 
Introducing the mean-field potential 

$(r) = / d^r' p{r')u{r-r'), (80) 



produced by the smooth density p{r) — APi(r), we obtain 



ap 5$ 

^(r) = -/3Pi(r) — (r). (81) 
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After integration, we find that the density of particles is given by the Boltzmann distri- 
bution 

p(r) = Ae-''*^, (82) 

where $ is related to p via Eq. (|80|l . Hence, Eq. H82I) must be regarded as an integro- 
diffcrential equation. 

We now consider the canonical distribution in phase space 

P^(ri,vi,...,r^,,v^) = -i-e-^{^-++^(---«)}. (83) 
Z(/3) 

Noting that the velocity distribution is Gaussian, we find that Eq. (|82(l is now replaced 
by the Maxwell-Boltzmann distribution 

/(r,v) = Ae-'^f'^+^W), (84) 

where, again, $ is related to p via Eq. 1)80(1 . The distribution function ((84(1 can also be 
obtained by minimizing the Boltzmann free energy FbIJ] at fixed mass M and tem- 
perature T (see Sec. 3). This method provides a condition of thermodynamical stability 
> 0, which is not captured by the equilibrium BBGKY hierarchy. Indeed, the ther- 
modynamical stability is related to the dynamical stability with respect to the iV-body 
Fokker-Planck equation (see Appendix B) corresponding to the time dependent BBGKY 
hierarchy (see also '15'). 

Appendix B: Brownian particles in interaction 

Bl. Stochastic processes in physical space In this Appendix, we derive the mean-field 
Fokker-Planck equation ((15(1 describing a gas of Brownian particles in interaction, using 
the method of Martzel & Aslangul tSQj- First, we consider a system of Brownian 
particles in physical space. The state of the system is completely specified by the A^-body 
distribution function P/v(ri, r^v, t) which gives the density probability of finding at 
time t the first particle in ri, the second in r2 etc. We assume that the evolution of the 
system is governed by the stochastic process 
dv 

— ^ = -/iV,C/(ri, ...,rAr) + y^M;{t), (85) 
at 

where Ri(i) is a white noise satisfying (Ri(i)) = and {Ra,i{t)Rbj{t')) — 6ijSabS{t — t'), 
where a, 6 = 1, D refer to the coordinates of space and i,j = 1, N to the particles. 
The particles interact via the potential [/(ri, ...,rjv) = J2i<j ~ ^j)- 
For a general Markovian process, we have 

Pw(ri,...,rw,i + At) = J d^(Ari)...d^(Arjv) Pw(ri-Ari,...,rjv-Arw,i) 

xu;(ri - Ari, ...,rN - ArAr|Ari, ArAr), (86) 

where w is the transition probability. Expanding the right hand side in Taylor series and 
introducing the Kramers-Moyal moments 

M„i...„„(ri,...,rjv) =limAt^o-r7 — 7 T / {Ari)...d^ {Atn) 

Ar ni\...nN\ J 

x(-Ari)"^..(-ArAr)""u;(ri,...,rjv|Ari,...,Arw), (87) 
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we get the N-body Fokker-Planck equation 



ni...niv 

where the sum runs over aU indices such that n-i > 1- For the stochastic process (jHSJ, 
only a few moments do not vanish, namely 

Mo...n,=i...o = MV»C/(ri,...,rAr), (89) 



Mo...„.=2...o = D,. (90) 
Substituting these results in the A^-body Fokker-Planck equation H88|l , we obtain 



dPi 



N 



dt 



N 

E 



_d_ 



D 



N 



dvi 



d 

^iPn^U{yi,...,yn) 



(91) 



BP /" B 

D,-^+ti{N-l) / d^r2P2(ri,r2,t)— u(ri-r2] 
dri J Byi 



(92) 



We note that the stationary solutions of this equation correspond to the canonical equi- 
librium distribution (|67|l provided that /i = (Einstein relation). Therefore, a gas 
of Brownian particles in interaction has a rigorous canonical structure. Integrating over 
r2,...,rjv, we obtain 
BPi _ 3 
Bt Br I 

Now, implementing a mean-field approximation which is valid in the proper thermody- 
namic limit TV > 1 with l3Nu ~ 1 jIHl: 

F2(ri,r2,i) =Pi(ri,t)Pi(r2,t) + 0(l/iV), (93) 
the foregoing equation can be rewritten 

BPi _ B 
Bt Bvi 



' BPi „ a* 
D,.— +fiPi — 
OTi Bri 



(94) 



which is clearly the same as 



Bp 
Bt 



V[D*(Vp + /3pV$)], 



(95) 



where $(r,t) is related to p{r,t) as in Eq. (|8U|I . 



B2. Stochastic processes in phase space We now consider the generalization of the 
preceding problem to the case of A^-coupled Brownian equations in phase space described 
by the stochastic process 



V,C/(ri,...,rAr) + \/2i?R,(i), 
where Ri(i) is a white noise. We start from the Markov process, 

Pjv(ri,vi,...,rAr,VA,,i + AO = j d^(Ari)d^(Avi)...d^(ArAr)d^(Avjv) 

xPAr(ri - Ari,vi - Avi, ...,rAr - ArAr,VAr - AvN,t) 
xw{ri - Ari,vi - Avi, ...,rN - Avn^^n - Avjv|Ari, Avi, ArAr, Avat). 



(96) 



(97) 
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Now, the transition probability can be written 

zz;(ri,vi, ...,rjv,VAr|Ari, Avi, Afat, Avat) 
(5(Ari - VI At)...(5(Arjv - VArAi)V'(ri, vi, tat, vw| Avi, Avjv). (98) 
The integration over Ari,...,Arjv is straightforward and yields 

Pw(ri,vi,..,rAr,vjv,i + Ai) = J d^(Avi)...d^(AvA,) 

xFAr(ri - viAt, vi - Avi, ...,rN - VArAt, vjv - Avjv,t) 
xipiri - viAt,vi - Avi, ...,rN - VNAt,VN - AvAr|Avi, Avat), (99) 
which is equivalent to 

PAf(ri+ViAt,Vi,...,rAr + VArAi,VAr,i + At) = J d^(Avi)...d^(AvAr) 

xPAf(ri, vi - Avi, ...,rjv, VAT - Avtv,*) 
x^(ri,vi - Avi, ...,r7v,v7v - AvAr|Avi, Avat). (100) 

Expanding the right hand side in Taylor series and introducing the Kramers-Moyal mo- 
ments 

M„^...„„(ri, vi, ...,rAr,VA') = limAt^o-n — T T / d^(Avi)...d^(AvAr) 

Ar ni\...npf\ J 

x(-Avi)"i...(-Avw)""V'(ri,Vi,..,r^,VAr|Avi,..,Av^), 
we obtain the N-body Fokker-Planck equation 



dP, 



N 



N 



dt 



dP, 



N 



dvi 



E 



Qn 



QriN 

-^{^^ni...nr,PN}, 



(101) 



(102) 



N 



i—1 " ni...njY 

where the sum runs over all indices such that J2i ni > I. For the stochastic process (|96|l. 
only a few moments do not vanish, namely 

Mo.. .0 = Cv, + V,C/(ri, r^), (103) 



M0...n,=2...0 = D. (104) 

Substituting these results in the A^-body Fokker-Planck equation p02|l , we get 



dp 



N 



dt 



N 

E 



dP, 



N 



dP, 



N 



Or, 



y— 



D^+^P^... 



(105) 



where — — ViC/(ri, tat). We note that the stationary solutions of this equation 
correspond to the canonical equilibrium distribution H83|) provided that ^ = Df3 (Einstein 
relation). Integrating over r2, V2, tat, vat, we obtain 

N 



dPi 



d 
9vi 



dPi 



(106) 
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Now, 

i=2 ^ j=2 i=2 ^ ^ 

^-{N-1) J d^r2d^V2 ^(ri-r2)|^(ri,vi,r2,V2,i). (107) 

Implementing a mean-field approximation which is valid in the proper thermodynamic 
limit iV > 1 with PNu ^ 1 48 : 

P2 (r 1 , VI , r2 , V2 , t) = Pi (r 1 , VI , t)Pi (r2 , V2 , i) + 0( 1/A^) , (108) 

we find that 



h Vi h (F)i = 

dt dri dvi dvi 



avi 



where 



(109) 

F)i = -nJ rf^r2rf^V2 ^(n -r2)Pi(r2,V2,t). (110) 



This can be rewritten 



df df ,^,df d 
— + V— + (F) — = — 
dt dr dv dv 



where 



(F) = -Va> = - 1 d^r' p{r', t)^{r- r'), 
is the mean-field force created by the particles. 



(Ill) 
(112) 
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